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We studied the interaction between a single hydrogen atom and a single graphene using classical molecular 
dynamics simulation with modified Brenner REBO potential. Three interactions, which are adsorption, reflec- 
tion, penetration, were observed. Overhang structure appears and creates an adsorption site on the backside of 
the graphene. It is considered that backside adsorption occurs under the two conditions that an incident hydrogen 
atom should have incident energy which is larger than the potential barrier of a hexagonal hole of the graphene 
and that after the hydrogen atom passes through the graphene, it does not keep its kinetic energy to be trapped 
by the adsorption site. The conditions explained that as the incident energy increased, the incident point of the 
backside adsorption shifted to the periphery of a hexagonal hole of the graphene in the simulation. Moreover, 
when a hexagonal hole of the graphene was expanded by the hydrogen atom incidence to the periphery of the 
hexagonal hole, its potential barrier was reduced. 
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1 Introduction 

In the context of research into nuclear fusion, the plasma 
surface interaction (PSI) problem has been studied [1-5]. 
A portion of the plasma confined in an experimental device 
falls onto a divertor wall, which is shielded by graphite 
or carbon fiber composite tiles. The incident hydro- 
gen plasma erodes these carbon tiles in a process called 
chemical sputtering. The erosion produces hydrocarbon 
molecules, such as CHx and C2Hx, which affect the plasma 
confinement. 

To solve the PSI problem, the mechanism of the 
graphite erosion has been researched using molecular dy- 
namics (MD) simulation [6-9]. Previously, we inves- 
tigated the PSI of graphite surface which consists of 8 
graphene sheets using the modified Brenner reactive em- 
pirical bond order (REBO) potential [10]. Subsequently, 
the isotope dependence of incident hydrogen atoms is in- 
vestigated [11]. These simulations achieve steady state 
of the graphite erosion, the incident energy linear depen- 
dence of total carbon yield accords with experimental re- 
sults [12,13]. 

The MD simulation of graphite surface showed that 
if incident energy was 5 eV, almost of all incident hydro- 
gen atoms were adsorbed by the graphite surface, while 
if the incident energy was 15 eV, most incident hydro- 
gen atoms were reflected. This adsorption and reflection 
can be explained by the MD simulation of the elemental 
processes which is the chemical reaction between a sin- 
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gle hydrogen atom and a single graphene [14-16]. This 
MD simulation indicates also that in a certain incident en- 
ergy, the hydrogen atom can be adsorbed to the backside 
of the graphene. On the other hand, the other MD simu- 
lation for graphene erosion due to hydrogen atom gas im- 
plied that para-overhang configuration which has hydro- 
gen atoms onto both side of the graphene is created and its 
C-C bond is easy broken [17]. Therefore, we consider that 
if the incident energy which brings about the backside ad- 
sorption is selected, the MD simulation of graphite surface 
shows a different dynamics. 

Here, we attach weight to not only the macroscopic 
simulation of graphite surface but also understanding a el- 
emental mechanism and dynamics. We, therefore, investi- 
gate the backside adsorption of the interaction between a 
single hydrogen atom and a single graphene in the present 
paper We describe the simulation model and method in 
S In ^ we present and discuss the simulation results. 
This paper concludes with a ^ 

2 Simulation Method 

A graphene [18] consists of 160 carbon atoms measuring 
2.13 nm x 1.97 nm, and is placed at the center of simula- 
tion box parallel to x-y plane. The simulation box in the 
X- and y-directions measures 2.13 nm x 1 .97 nm with pe- 
riodic boundary condition. In the z-direction, we do not 
prepare the boundary of the simulation box and the initial 
coordinate of the center of mass of the graphene is z = A. 
The initial graphene temperature is set to K. The carbon 
atoms of the graphene are relaxed initially and compose 
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and backside adsorption where the graphene temperature 
is OK. 



Fig. 1 The incident energy dependence of the rates of adsorp- 
tion, reflection and penetration where the graphene tem- 
perature is K. Dash-dotted line with open circle, long- 
dashed line with filled triangle, and short-dashed line 
with square denote the rates of adsorption, reflection and 
penetration, respectively. 



completely flat graphene structure. A hydrogen atom is 
injected from z = 4 A parallel to the x-y plane, that is, nor- 
mal incidence to the graphene surface. Incident position in 
the X- and y- directions are determined under a uniformed 
distribution function. The incident energy Ei decides an 
initial momentum vector (0, 0, pu) as follows: 



-\/2mEi. 



(1) 



where and m is the mass of the hydrogen atom. 

We performed our MD simulation under NVE condi- 
tions, where the number of atoms, volume, and total energy 
are conserved. The simulation time was developed using 
second order symplectic integration [19]. The chemical in- 
teraction was represented by the modified Brenner REBO 
potential [15,20]: 



i,j>i 



■bij{{r]Ae^}Ae^''})vL(nj) 



(2) 



where r,j is the distance between the /-th and j-th atoms. 
The functions V^j^ and V^j^ represent repulsion and attrac- 
tion, respectively. The function bij generates multi-body 
force. To conserve the accuracy of the calculation, the time 
step was 5 x 10"'** s. 

We note the condition to judge the type of reaction, 
which corresponds to the condition to finish the simula- 
tion. While the incident hydrogen atom interacts with one 
carbon atom, we count a trapped time. When the hydrogen 
atom leaves the carbon atom or starts interaction with the 
other carbon atom, the trapped time is cleared. When the 
trapped time reaches 0. 1 ps (20000 time steps), the simula- 
tion is finished and this reaction is regarded as adsorption. 



Moreover, if the relative position from the nearest carbon 
atom to the hydrogen atom in the z-coordinate is positive, 
the reaction is front adsorption, while if the relative posi- 
tion in the z-coordinate is negative, the reaction is backside 
adsorption. When the hydrogen atom leaves the nearest 
carbon atom, we begin counting a escaping time. If the 
hydrogen atom interacts a carbon atom again, we clear the 
escaping time and restart counting the trapped time. When 
the escaping time reaches 0.05 ps (10000 time steps), the 
simulation is finished. And, if the escaping hydrogen atom 
has positive momentum in the z-coordinate, this reaction 
is reflection, while it has negative momentum, the reaction 
is penetration. 

We repeated the above simulation 200 times for each 
incident energy, where the incident position in the x- and 
y- directions are every time changed randomly under the 
uniformed distribution function. Since we counted the 
types of the reactions, we obtained the rates of the reac- 
tions. Of course, the sum of the rates of the adsorption, 
reflection and penetration is always 1 . 

3 Results and Discussion 

The simulations were executed for incident energy of 0. 1 
eV to 200 eV. The three interactions, which are adsorption, 
reflection, penetration, were observed. Figure[T]shows that 
the rates of the three interactions depend on the incident 
energy. Especially, the rate of adsorption has two peaks. 
Next, Fig. |2] shows the rates of the front and backside 
adsorption. Although the front adsorption occurs around 
the both adsorption peaks, the backside adsorption occurs 
around the high incident energy adsorption peak only. In 
the process of front adsorption, the nearest carbon atom 
which is connected to the incident hydrogen atom by a co- 
valent bond is pulled out from flat surface of the graphene. 
This is called a overhang structure and creates an adsorp- 
tion site for a hydrogen atom [14, 15]. In the backside ad- 
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The potential barrier of a hexagonal hole of the graphene. 
The mapped value is potential energy when the hydrogen 
atom is located in the plane of z = A. 



sorption also, the overhang structure appears. That is, an 
adsorption site is created in the backside of graphene after 
the incident hydrogen atom goes through the graphene. 

The incident energy dependence of the rates of inter- 
actions Fig. [T|is understood using energetics, which relates 
;r-electron on the graphene surface, structure transforma- 
tion to the overhang structure and the potential barrier of 
a hexagonal hole of the graphene [15]. In the present pa- 
per, we discuss the backside adsorption. We consider that 
the backside adsorption needs the following two condi- 
tions. One is that to go through the graphene, the incident 
hydrogen atom should have the incident energy which is 
larger than the potential barrier of a hexagonal hole of the 
graphene. The other is that after the hydrogen atom passes 
through the graphene, it is trapped by the adsorption site 
due to the overhang structure. Namely, to occur the back- 
side adsorption, the kinetic energy of the hydrogen atom 
is diffused to carbon atoms of the graphene and must not 
remain. Concerning the first condition, the potential bar- 
rier of a hexagonal hole of the graphene, where the hydro- 
gen atom is located at z = A, is shown by Fig. [3] The 
center of this figure corresponds to the center of a hexag- 
onal hole of the graphene. The six white regions indicate 
higher potential baiTier of more than 50 eV because car- 
bon atoms exist there. The point of minimum potential 
barrier is not the center of a hexagonal hole of graphene 
and appears between the locations of carbon atoms and 
the center of a hexagonal hole of the graphene. Figure |4] 
shows incident points and the types of their interactions. It 
is understood and consistent with the potential barrier of 
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Fig. 4 Incident points v.s. the types of its interactions. Circles 
indicate incident points which bring about hydrogen atom 
adsorption where the opaque circles correspond to back- 
side adsorption. Triangles and squares are reflection and 
penetration, respectively. Big six circles and six line rep- 
resent six carbon atoms and six covalent bonds compos- 
ing a hexagonal hole of the graphene. 



a hexagonal hole of the graphene Fig. [3] that the backside 
adsorption occurs around the points of minimum potential 
barrier and does not occur around the center of a hexago- 
nal hole of the graphene. As the incident energy increases, 
the incident points for the backside adsorption shift from 
the points of minimum potential baiTier to the periphery 
of a hexagonal hole of the graphene which is over the co- 
valent bonds between carbon atoms. This phenomenon is 
explained as follows. In the incident energy of 20 eV, be- 
cause all incidences become the front adsorption or reflec- 
tion, the hydrogen atom does not have enough incident en- 
ergy to go over the potential barrier of a hexagonal hole of 
the graphene. In the incident energy of 25 eV, because the 
incident energy is higher than but close to the minimum po- 
tential baiTier, only the hydrogen atoms which are injected 
around the points of the minimum potential barrier go 
through the hexagonal hole of the graphene and are trapped 
by adsorption site of the backside of the graphene. When 
the incident energy increases (27 and 30 eV), the region 
in which the potential barrier of a hexagonal hole of the 
graphene is lower than the incident energy extends. There- 
fore, the backside adsorption occurs in the periphery of a 
hexagonal hole of the graphene. However, the backside ad- 
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sorption around the potential minimum point changes into 
penetration because the hydrogen atom leaves the kinetic 
energy enough to escape out of the adsorption site after it 
passed through a hexagonal hole of the graphene. More 
important thing is that while the backside adsorption hap- 
pens in the periphery of a hexagonal hole of the graphene, 
it does not happen in the center of a hole of the graphene 
in spite of same potential barrier. It is demonstrated by the 
trajectory of atoms in the MD simulation that when the hy- 
drogen atom is injected into the periphery of a hexagonal 
hole of the graphene, near two carbon atoms move and the 
hexagonal hole of the graphene is expanded. The expan- 
sion of the hexagonal hole of the graphene seems to reduce 
the potential barrier of the hexagonal hole of the graphene 
because the distance between hydrogen atom and the near 
carbon atoms become larger. On the other hand, if the hy- 
drogen atom is injected to the center of a hexagonal hole of 
the graphene, it must push six carbon atoms to expand the 
hexagonal hole of the graphene. It is considered that at the 
incident energy of 30 eV, even if the hydrogen atom can 
move two carbon atoms, it cannot move six carbon atoms 
because of their larger mass. Consequently, the backside 
adsorption and penetration hardly occur in the center of a 
hexagonal hole of the graphene. 

4 Summary 

We study the interaction between a single hydrogen atom 
and a single graphene using classical molecular dynamics 
simulation with modified Brenner REBO potential. The 
three interactions, which are adsorption, reflection, pene- 
tration, were observed. Especially, in the present paper, 
we discuss the backside adsorption. As well as the front 
adsorption, the overhang structure appears in the backside 
of the graphene and creates the adsorption site. It is consid- 
ered that the backside adsorption occurs under the two con- 
ditions that the incident hydrogen atoms should have the 
incident energy which is larger than the potential barrier of 
a hexagonal hole of the graphene and that after the hydro- 
gen atom passes through the graphene, it does not keep its 
kinetic energy to be trapped by the adsorption site. The hy- 
drogen atom in the incident energy of 25 eV brings about 
the backside adsorption only around the points of the mini- 
mum potential barrier, which is not the center of a hexago- 
nal hole of the graphene and is the region between the cen- 
ter of a hexagonal hole and the locations of carbon atoms. 
When the incident energy increases, the incident points for 
the backside adsorption shift to the periphery of a hexago- 
nal hole of the graphene. In addition, the backside adsorp- 
tion around the potential minimum point change into pen- 
etration because the hydrogen atom leaves enough kinetic 
energy to escape from the adsorption site after it passed 
through a hexagonal hole of the graphene. Moreover, when 
a hexagonal hole of the graphene is expanded by the hydro- 
gen atom incidence to the periphery of the hexagonal hole. 



its potential barrier is reduced. However, even if the hy- 
drogen atom can move two carbon atoms in the incidence 
to the periphery of the hexagonal hole, it cannot move six 
carbon atoms in the incidence to the center of the hexago- 
nal hole. Therefore, the backside adsorption and penetra- 
tion hardly occur in the center of a hexagonal hole of the 
graphene. 
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